Metal-insulator transition in hydrogenated graphene as 
manifestation of quasiparticle spectrum rearrangement of 

anomalous type 

Yuriy V. Skrypnyk 1 and Vadim M. Loktev 2 

1 G. V. Kurdyumov Institute of Metal Physics, 

National Academy of Sciences of Ukraine, 

Vernadsky Ave. 36, Kyiv 03680, Ukraine 
2 Bogolyubov Institute for Theoretical Physics, 

National Academy of Sciences of Ukraine, 
Metrolohichna Str. 14-b, Kyiv 03680, Ukraine 

Abstract 

We demonstrate that the spectrum rearrangement can be considered as a precursor and a basis for 
the metal-insulator transition observed in graphene dosed with hydrogen atoms. The Anderson- 
type transition is attributed to the Fermi level's entering into the quasigap, which develop in 
the vicinity of the impurity resonance energy due to the anomalous spectrum rearrangement. 
Theoretical results for the Lifshitz impurity model are in a reasonable agreement with available 
experimental data. 
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I. INTRODUCTION 



Since the existence of Dirac quasiparticles has been proved for graphene, one of the most 
intriguing issues of its physics is the possibility of their localization by whichever imperfection 
that can appear in the honeycomb lattice." Moreover, It should not be overlooked that the 
effect of disorder on the transport of quasiparticles is sensitive to the nature of inhomogeneity, 
be it caused by short-range or long-range defects, neutral or charged impurities, adsorbed 
or interstitial atoms or molecules, vacancies, spatial distortions or structural irregularities, 
including ripples on the graphene sheet and other long-wave random modulations. All these 
different types of imperfections necessitate corresponding dedicated studies, which can not be 
accomplished by resorting to a single impurity model. Since these imperfections are rather 
dissimilar both in their character and action on basic properties of graphene, respective 
theoretical models should be diverse as well. 

Early experiments on graphene-based devices, which were engineered like a commonplace 
field effect transistors, revealed that the sample conductivity never drops below a certain 
minimal value.- This fact, indeed, considerably reduced audacious expectations that corre- 
sponding devices are capable to serve as next-generation electronic switches. Nevertheless 
(or, maybe, exactly because of this), the minimal conductivity existence has produced quite 
a stir, and its origin has been relentlessly debated. Distinctive features of charge carriers, 
which were shown to obey the linear dispersion, constituted the core of this discussion. 
Uniqueness of the electron subsystem in graphene were pushed to its limits so much that 
former physics of semiconductors were sometimes categorically declared being utterly un- 
suitable for this uncommon material. It has been speculated that massless, according to 
their Dirac dispersion, charge carriers can not be localized by any degree of disorder caused 
by lattice imperfections or impurity centers. The presumed impossibility to localize Dirac 
excitations were directly linked to the minimal conductivity phenomenon by a simple reason- 
ing. Since the product of the wave vector modulus and the mean free path is confined from 
below for propagating states according to the Ioffe-Regel criterion, the conventional Drude 
expression immediately yields the minimum value of ~ e 2 /h for each conducting channel. 
This argument, along with available experimental results, has led to an opinion that the 
minimum conductivity value has a universal character for graphene, and is expressed only 
through fundamental constants. 
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In the mean time, there appeared occasional theoretical studies, which did not deny a 
possibility to localize charge carriers in graphene, and argued for the mobility edge appear- 
ance under certain circumstances, either for a disorder of a general type,— >2 or for a specific 
impurity model, in particular, for substitutional defects,— and for the Anderson model of 
disorder.—"— Later, it has become apparent that the minimal conductivity in graphene does 
depend on the sample quality and noticeably varies with the impurity concentration.—"— 
Furthermore, not so long ago, a metal-insulator transition (MIT) in graphene dosed with 
atomic, which is essential, hydrogen has been convincingly observed. The MIT has been 
manifested by an increase of the room temperature resistance in about four orders of 
magnitude.— The transition has been reported for the grown on SiC surface graphene with 
a low hydrogen coverage, namely, around 0.1%. The presence of the mobility gap has been 
also reported for the fluorinated and ion bombarded graphene.— >2i 

As well, corresponding data of ARPES measurements indicated the disappearance of 
quasiparticle excitations in the system. The MIT of the Anderson type comes about when 
the Fermi level finds itself inside a domain of localized states. Since the graphene in the ac- 
tual experiment was pristine before the hydrogen deposition and possessed sounding metal- 
lic properties with conventional Fermi-liquid behavior, it is tempting to anticipate that 
some amount of deposited hydrogen is sufficient to open up a reasonably wide region filled 
with localized states, or a quasigap, in the spectrum. In essence, this draft description of 
the leading to the MIT process corresponds to the well-known phenomenon of spectrum 
rearrangement , 281 which deals with decisive modifications in the elementary excitation 
spectrum upon increasing the impurity concentration in the system. When the quasigap 
is located in the vicinity of the Fermi level position, the MIT can take place in spite of 
the rather low hydrogen concentration. Below we are going to link together the predicted 
spectrum rearrangement in graphene with point defects and the observed MIT under the hy- 
drogen dosing. 19 The spectrum rearrangement takes place when a single impurity induces a 
local or a resonance level in the spectrum. It has been demonstrated recently that graphene 
almost inevitably contains traces of impurity centers, which are capable in producing reso- 
nance states.^ 
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II. CRITIC ALITY AND TYPES OF SPECTRUM REARRANGEMENT 



Let us discuss some general spectral properties of impure systems and assume for this 
purpose that inside an isolated host band in a hypothetical d- dimensional system a single 
impurity center accounts for a resonance state with the energy e r > 0, which is measured 
against the band edge and has been made dimensionless by the corresponding bandwidth. 
The dispersion relation for this host band is supposed to be of the form e(k) ~ k p . Then, 
the spatial behavior of the host Green's function at a given energy inside the band, apart 
from any other details, should be governed by oscillations with the characteristic radius 
r ~ £~ 1 / p , which is expressed in units of the appropriate lattice constant. At the energy e r , 
these oscillations determine the effective radius r im ~ e^ 1 ^ of the impurity state. Through 
a closeness of the resonance energy to the band edge, the effective radius of the impurity 
state may far exceed the lattice period. As the impurity concentration is increased, the 
individual impurity states are becoming more tightly packed. It can be anticipated that a 
significant spatial overlap between these states is capable in provoking qualitative changes in 
the spectrum. The average distance between impurities f ~ c -1 ^, where c is the (relative) 
impurity concentration, is gradually shrinking with increasing c, and at some stage of the 
process becomes of the order of rj m . This general condition defines the critical concentration 
csr ~ £r of the spectrum rearrangement in the corresponding impure system. It seems 
justifiable to reiterate here that the value of csr can be fairly low due to the smallness of 
the energy e r . 

Graphene is, evidently, a two-dimensional system, so that d = 2, and features the linear 
dispersion of charge carriers, which implies that p — 1. Consequently, the critical concen- 
tration of spectrum rearrangement in the impure graphene is expected to be ~ ef, where 
energy is counted from the Dirac point, at which the valence band and the conduction band 
coincide. In a case of the conventional dispersion with p = 2, the identical in appearance 
relation between e r and csr is inherent in four-dimensional systems. In this sense graphene, 
as regards the spectrum rearrangement, can be formally viewed as a conventional system of 
increased spatial dimensionality. This fact further accentuates the uniqueness of graphene, 
since the spectrum rearrangement has not been studied for such systems. 



Let us return back to experimental data contained in Ref. Il9|. The comparison of angle 



integrated spectrum for clean and hydrogen dosed graphene suggests that the dimensional 
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resonance energy E r lies somewhere around 0.25 eV lover than the initial Fermi level, while 
the Dirac point position Ed in the undoped crystal is about 0.45 eV below it. A bit more 
careful analysis of the spectrum rearrangement in impure graphene,— which is based on the 
expansion of the self-energy into a cluster series,— brings about the following expression for 
csr, which, for a convenience, is rewritten here in terms of the hydrogen coverage, 

(E r — Ed) 2 n , 
n SR ~ n c — , (1) 

where 

4 

n c = ~ 3.8 x 10 15 cm~ 2 (2) 

V3a 2 

is the areal density of carbon atoms in graphene, a = 0.246 nm - the lattice constant 
of graphene, ngR - the critical coverage of hydrogen atoms, W ~ 6.3 eV - the bandwidth 
parameter, and the logarithmic correction is omitted because of the relatively large resonance 
energy. Substitution of the guessed value for the resonance energy results in usr that is 
reasonably close to the reported MIT critical hydrogen coverage of 3.8 x 10 12 cm -2 .— As 
a matter of fact, the critical coverage of the spectrum rearrangement nsR, should not be 
identical to the experimentally obtained hydrogen coverage required for the MIT, which will 
be discussed later. 

It should be also noted that there are two main types of the spectrum rearrangement: 
the cross one, which is usually accompanied by a sharp single impurity resonance, and the 
anomalous one. The first type of the spectrum rearrangement results in a quasiparticle dis- 
persion that looks like a familiar hybridization between the host branch and a dispersionless 
branch at the impurity resonance energy. In this case, two new branches are separated by a 
gap, which widens with an increase in the impurity concentration. Consequently, two differ- 
ent wave vectors correspond to the same energy in this split spectrum. However, this feature 
has not been detected in the ARPES measurements on the hydrogenated graphene. The 
second type of the spectrum rearrangement, which is usually encountered in low- dimensional 
systems, is of a more diffused nature and frequently corresponds to a considerably smeared 
single impurity resonance. This anomalous spectrum rearrangement is characterized by 
the opening of the quasigap, which is filled with localized states and separates two non- 
overlapping branches of extended states exhibiting a renormalized dispersion. 
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III. INTERPLAY BETWEEN SPECTRUM REARRANGEMENT AND ANDER- 
SON TRANSITION 



The anomalous type of the spectrum rearrangement has been shown to unfold in graphene 
for impurities described by the Lifshitz model. 25 Within this model, identical impurity cen- 
ters are supposed to be randomly distributed on sites of the host lattice. Each impurity is 
only allowed to change the on-site energy at its location in the corresponding tight-binding 
Hamiltonian. It is well-known that the electron subsystem in graphene encompasses prac- 
tically free 7r-electrons, which number is equal to the number of lattice sites. As soon as an 
additional hydrogen atom lands on graphene, its uncoupled s-electron immediately enters 
a chemical bond with one of the 7r-electrons. That should make the latter one localized 
almost completely at the impurity site. In the first approximation, one can assume that the 
resulting perturbation comes only from the appearance of strong attracting potential on the 
occupied carbon atom. Thus, the disordered system is described by following operators, 

11 11 H im , H im = V L ^2 w4a c «a> (3) 

na 

where 

H = t c Lc m/3 (4) 

na,m/3 

is the host Hamiltonian, in which summation is restricted to nearest neighbors, n runs 
over lattice cells, indices a and (3 enumerate sublattices, c^ na and c na are the creation and 
annihilation operators at the corresponding lattice site, t ~ 2.7 eV - is the tight-binding 
transfer integral for the 7r bands in graphene, the variable rj nQ takes the value of 1 with 
the probability n H /ric or the value of with the probability 1 — n H /ric, where tih is the 
hydrogen coverage, and Vl < is the difference in the on-site potential at the defect position. 

Within this model, the local density of states (LDOS) at the lattice site occupied by a 
hydrogen reads, 

where go{E) is the diagonal element of the host Green's function 

g = (E-H )-\ (6) 
In the vicinity of the Dirac point in the spectrum, namely, within the window stretching up 
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to around 0.75 eV to each side of it, this diagonal element can be easily approximated, 



E — Ed ( E — Ej)\ . \E — Ed\ 



gJE)&- — =^ln - — — — i7r— — — (7) 

v ; w 2 \ w I w 2 



where 



W = A/ nV3t (8) 



is the same bandwidth parameter as in ([I]). Substitution of this approximation to Eq. (jSJ) 
yields 

\£\ 

PUE) W W[(l-2^1n|^ + ( 7 r^)T (9) 

where 

S — (E — E D )/W. (10) 

The LDOS at the impurity site is shown in Fig. [T] for Vj, = —25 eV. As follows from the 
Figure, the resonance energy E r is located approximately 0.2 eV above the Dirac point and 
thus corresponds to the experimentally observed peak. While the required change in the 
on-site potential is quite substantial, the large perturbation magnitude should be considered, 
first of all, as a result of the chosen simplified impurity model. Unless another is explicitly 
specified, the perturbation value of —25 eV will be used below for all subsequent estimations. 

When a certain amount of deposited hydrogen atoms is taken into account, the negative 
potentials induced by them on the occupied lattice sites impede the electron movement. For 
the Lifshitz impurity model, which has been described above, the course of the spectrum 
rearrangement has been studied both analytically and numerically.— 1 ^ The analysis has 
shown that states are subject to localization near the resonance energy, where the impurity 
scattering is the strongest. At some coverage, a quasigap occupied by localized states opens 
around the energy E r . With increasing impurity concentration, the quasigap gradually 
broadens and the degree of localization inside it rises. This quasigap expansion and the 
localization enhancement are not symmetric about the energy E r and are more expressed 
above the resonance (for an attractive potential). Finally, the quasigap consumes all the 
space between E r and Ejj in the spectrum, which corresponds to the onset of the spectrum 
rearrangement by the definition. In other words, the spectrum rearrangement does follow 
the anomalous scenario. 

In the experiment discussed, dosed hydrogen atoms act not only as scattering centers, 
but also as acceptors. Thus, their presence inevitably lowers the Fermi level of the system. 
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FIG. 1: Local density of states at the impurity site for Vl = —25 eV. 

Therefore, while the hydrogen coverage is increasing, the Fermi level and the upper boundary 
of the quasigap, where the localization of states is most pronounced, are moving toward each 
other. Sooner or later, the Fermi level should appear inside the quasigap, which will cause 
the MIT of the Anderson type. Since the energy E r is positioned somewhere in-between E D 
and the bare Fermi level Ep are \ the critical coverage of the spectrum rearrangement should 
be close but not identical to the critical coverage of the MIT. 

The acceptor effect of the deposited hydrogen atoms is also implicitly contained in the 
Lifshitz model. Strong attracting potential Vl < produces a deep local level below the 
conduction band. Because of its remoteness from the continuous spectrum, the correspond- 
ing narrow impurity band should hold almost UhN c / '(2nc) energy levels, where N c is the 
total amount of carbon atoms in the system. Thus, taking into account the spin degeneracy, 
this impurity band consumes nearly one electron per deposited hydrogen atom. However, 
apart from these qualitative considerations, the adopted impurity model allows to address 
the question on the Fermi level position more precisely. Consider the (normalized) total 
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number of states with energies that are less than a specified one, 

N(E)= [ E p(E)dE, (11) 



where p(E) is the density of states. In a disordered system this quantity can be expanded 
in powers of the impurity concentration,— 

N(E) = N (E) + — ^ arg(l - V L g (E)) + . . . , (12) 
ixn c 

where the presence of two sublattices in graphene and the Lifshitz impurity model are taken 
into account, and N (E) is the total number of states in the host system. Keeping in mind 
that the deep impurity band is modeling the acceptor effect of impurities, it is natural to 
demand the conservation of the number of occupied states in the system with varying the 
amount of dosed hydrogen. This yields a kind of a balance condition, 

E F — E D \ (E F — E D ) 2 2n H (Ep are ^ — E D ) 2 



sgn ) Z^=*- + ar g(1 - Vm(EF)) = ~ ^ , (13 ) 

where the constant term is dropped from the both sides of the equation, and Eq. (j7J) is used 
to obtain Nq(E) for the clean graphene. From here and on this condition will be employed 
to determine the Fermi level position at a given dosing level. It is worth mentioning that 
the Dirac point position (which should not be confused with its value for the unperturbed 
system Ed) is gradually shifted with increasing the impurity concentration by approximately 
^//V^/^c- Since the Dirac point and the Fermi level are both moving in one direction, namely 
to lower energies, the distance between them does not shorten remarkably. 



IV. CONDUCTION BAND SPECTRUM 



For njj < nsRi the presence of attracting impurities does not distort significantly the 
quasiparticle dispersion in the valence band. In contrast, impurities noticeably affect the 
conduction band. The self-energy X, which enters the Dyson equation 

G = g + gHG (14) 

for the averaged over impurity distributions Green's function of the disordered system 

G = {(E-H)- 1 ) (15) 
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can not be attributed only to the impurity scattering in the actual experiment of Ref. [19, 
because the linewidth of the ARPES data is not negligible even at the absence of impurities. 
To deal with this issue, we assume that all other sources of its broadening, except the 
controlled amount of hydrogen atoms, can be roughly described by a constant dumping 
term, 

£«£ iTB -iri, r>o, (16) 

where Sj m is those part of the self-energy, which results from the impurity scattering, and 
J is the identity matrix. A comparison of the density of states obtained by the numerical 
simulation with the one resulted from the standard average T-matrix method showed that at 
n H < n SR this approximation works rather well in the vicinity of the impurity resonance.— 
Within this single-site approximation, the impurity-induced part of the self-energy is diag- 
onal in lattice sites and sublattices for point defects, 



an 



nc I. (17) 



1 - (1 - ^)V L g (E) 



Thus, taking into account both assumptions, the total self-energy is also diagonal, 

£ « E(E)I. (18) 
According to the conventional expression for the spectral function, 

A(E, k) « -- Im , (19) 



where 



is the Fermi velocity, the wave vector k is counted from the corresponding Dirac point, and 
only one branch, which is manifested in the actual experiment, is retained, it is straight- 
forward to put down corresponding expressions for the inverse width of the momentum 
distribution curve at its half-height, 

L ^ ~ "SEW <21) 



and for the inverse of its median 



- E^m- (22) 
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FIG. 2: The mean free path Lp (solid line) and the inverse wave vector k F l (dashed line) at the 
Fermi level vs the hydrogen coverage. 

Both this quantities are taken at the Fermi level, Lp = L(Ep) and hp 1 = l/k(Ep), and 
depicted in Fig. |2] in their dependence on the hydrogen coverage. The Fermi level position 
is calculated from Eq. (fT3|) . and the external damping T is chosen to be 0.03 eV. 

The ARPES linewidth analysis given in Fig. 3(c) of Ref. |l9| can be compared with Fig. |2j 
Since the inverse momentum width L(E) corresponds to the magnitude of the mean free 
path of charge carriers, the Ioffe-Regel criterion^ is violated for states at the Fermi level, 
i.e. kpLp = 1, at somewhat higher hydrogen coverage, namely, around 9.5 x 10 12 cm -2 , 
than in the experiment - 6.5 x 10 12 cm~ 2 . This discrepancy, for the most part, arise from 
the noticeably reduced wave vectors in the conduction band of the sample containing no 
hydrogen. The reduction is evident in comparison with the dispersion in the valence band, 
or with the idealized graphene fll]). The experimentally obtained quasiparticle dispersion 
looks already distinctly distorted before any hydrogen deposition. The distortion can be 
attributed to the possible presence of unidentified impurities in the system.— To a certain 
extent, this issue can be crudely addressed by a corresponding adjustment of the hopping 
integral t in the Hamiltonian. Results of such quick guesstimate are shown in Fig. [3j 
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FIG. 3: The mean free path Lp (solid line) and the inverse wave vector kp 1 (dashed line) at the 
Fermi level vs the hydrogen coverage for t = 3.5 eV, Vl = —35 eV, and T = 0.04 eV. 

The deliberate adjustment of the transfer integral t pushes kp 1 to where it should be at 
the zero hydrogen coverage. Correspondingly, the Ioffe-Regel criterion is violated at lower 
hydrogen coverage matching the experimental data. However, such a forthright measure 
will be suitable only in the Fermi level vicinity. It is more consistent to work out a proper 
model of the host system that is capable in simulating its experimentally obtained spectrum. 
Anyhow, the resulting adjustment of the zero coverage Fermi wave vector moves the entire 
curve for the inverse Fermi vector upwards and, thus, dimmish the impurity concentration, 
at which the Ioffe-Regel criterion ceases to hold. 

It should be mentioned that, according to the routine of the spectrum rearrangement 
analysis, one should expect the Fermi level to enter inside a domain of localized states at 
kpLp ~ v3/2. This will happen for the idealized host spectrum at the hydrogen coverage 
around 10.5 x 10 12 cm -2 , which is a little bit higher than the coverage corresponding to 
kpLp = 1. With further increase in the impurity concentration, the average T-matrix ap- 
proximation (1171) becomes non- valid because of the increased scattering on impurity clusters, 
and the approach outlined above is only suitable for signaling a strong localization at the 



12 




4 6 

n H xl0 12 , cm" 2 



8 



FIG. 4: Logarithm of the dimensionless zero temperature conductivity vs the hydrogen coverage. 

Fermi level. On the other side, the experimentally detected sharp increase in the sample 
resistance, see Fig. 3(b) of the Ref. lj], also succeeds the Ioffe-Regel criterion violation. 
Within the Kubo approach the zero temperature conductivity can be written as follows, 28 



&cond 



@condi 



&cond 



2 

7T 



1 + (cot if F + taXL(p F )(J^ - Lfp^j 



(23) 
(24) 



where 

ifp = Bxg(Ep — Tj(Ef)). (25) 

The dependence of the dimensionless conductivity a con d on the hydrogen coverage is plotted 
in Fig. HI It is evident from Fig. H]that the conductivity is falling down nearly exponentially 
with increasing the impurity concentration. Indeed, the Kubo formula and the average T- 
matrix approximation are becoming not so reliable as approaching the mobility edge. Thus, 
the zero temperature conductivity can not be satisfactorily described by Eqs. f[T7|) and (|2~4"|) 
close to the MIT. 

The concentration dependence of the energy distribution curve at the Fermi wave vector 
k F can be easily re-created for the Lifshitz impurity model with the help of Eqs. ([TBI . ffTBT) . 
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FIG. 5: Evolution of the energy distribution curve at Jzf with increasing the hydrogen coverage. 

(TT7|) . and (1T91 . The corresponding magnitude, A(E,kir), is given in Fig. as a function 
of the binding energy E — Ep and the hydrogen coverage Uh- It is clearly seen from the 
Figure that the energy width of the quasiparticle peak at /jp is gradually increasing with 
increasing the hydrogen coverage. Close to the critical impurity concentration for the MIT, 
the width of the energy distribution curve occupies nearly all conduction band below the 
Fermi level, which indicates the presence of a developed quasigap in the spectrum. This 
broadening translates into the complete breakdown of the quasiparticle picture near the 
MIT. The characteristic curve shape consisting of two competing peaks reproduces well the 
experimentally observable one, which can be seen in Fig. 2(c) of Ref. [l9|. At that, the peak, 



which is developing with increasing the hydrogen coverage, is connected with the impurity 
resonance. 

Finally, we render the density plots of the spectral function ( 1T9|) at hydrogen coverages 
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that are close to the critical one for the MIT in order to reproduce Fig. 1(d) of Ref. 
In Fig. [6] the Fermi level is about to enter the spectral region near the resonance energy, 
where the quasigap is forming due to strong impurity scattering. In Fig. [7] the Fermi level is 
already inside the developed quasigap. The pattern of Fig. [7] suggests the apparent upturn 
of the dispersion as approaching the Fermi level of the system, which has been noticed in 
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FIG. 6: Density plot of the spectral function at nu = 5 x 10 12 cm' 2 . 

the experiment. It is worth mentioning here that point defects does not cause a uniform 
broadening of the spectral function, on the contrary, the impurity-induced broadening is 
most pronounced in the vicinity of the resonance energy E r in this case. 

V. CONCLUSION 

To summarize, we have established that point defects are introducing a new length pa- 
rameter into the impure graphene. This length parameter results from the presence of a 
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FIG. 7: Density plot of the spectral function at nu = 10 x 10 cm . 

single-impurity resonance and shows up as the effective radius of a single-impurity state. 
When the average distance between impurities decreases up to this effective radius, the 
quasiparticle spectrum undergoes the cardinal rearrangement. The spectrum rearrange- 
ment is manifested in the opening of the quasigap around the impurity resonance energy 
in the spectrum. The quasigap progressively widens with increasing the impurity concen- 
tration. The Fermi level moves due to the doping effect of impurities, as well as confining 
the quasigap mobility edges are moving due to its expansion. If the Fermi level crosses one 
of the mobility edges and enters the quasigap, where states are localized, the MIT of the 
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Anderson type takes place in the disordered system. 

We have presented arguments, which confirm that the above scenario can be considered 
as a sounding candidate for the explanation of the experimentally observed MIT in graphen 
dosed with hydrogen. Even employing the simple Lifshitz model for the impurity centers, it 
appears possible to achieve a semi- quantitative interpretation of the experimental data. The 
MIT of the Anderson type in this case is prompted by the gradual lowering of the Fermi level 
due to the acceptor effect of the hydrogen atoms combined with the persistent raising of the 
mobility edge due to the ongoing spectrum rearrangement. Thus, the spectrum rearrange- 
ment acts as the main cause of the MIT, and, as a common phenomenon in semiconductors, 
provides the basis for understanding the physics of the process. Indeed, the oversimplified 
impurity model is not able to convey all the detail of the system's behavior. To fulfill this 
task, more sophisticated impurity models are required. Among them the two-parametric s-d 
model looks like a more natural choice for the deposited hydrogen atom. However, more so- 
phisticated impurity models should not change the general physics of the transition process, 
which is already captured by the Lifshitz impurity model. 
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